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We derive the photocount statistics of the radiation emit- 
ted from a chaotic laser resonator in the regime of single-mode 
lasing. Random spatial variations of the resonator eigenfunc- 
tions lead to strong mode-to-mode fluctuations of the laser 
emission. The distribution of the mean photocount over an 
ensemble of modes changes qualitatively at the lasing transi- 
tion, and displays up to three peaks above the lasing thresh- 
old. 

PACS numbers: 42.50.Ar, 05.45.-a, 42.25. Dd, 42.55.-f 



Nonlinear interactions in active optical media can dras- 
tically affect the statistics of the photon field. The text- 
book example is the single-mode laser Below the 
laser threshold, losses outweigh the gain by linear am- 
plification such that nonlinear saturation is negligible. 
Then the photon statistics is well represented by a ther- 
mal distribution, and the probability P„ for the occupa- 
tion of the mode with n photons decays as the power law 
Pn ~ [n/(l-l-fi)]" (n is the mean photon number). Above 
threshold nonlinear interactions stabilize the field inten- 
sity. The (relative) intensity fluctuations are strongly 
reduced below the value found for a thermal distribu- 
tion, and far above threshold P„ approaches the Pois- 
sonian distribution P„ = [n"/?^!] exp(— n) characteristic 
for a coherent state. These facts have been known since 
the development of the quantum theory of lasing in the 
1960s. 

Interest in the photon statistics of amplifying media 
was recently renewed by experiments on artificially fab- 
ricated random media. In random media radiation is 
scattered in an irregular, chaotic way. Practical realiza- 
tions include laser dye solutions |^ and semiconductors 
Q with randomly fluctuating dielectric constant as well 
as chaotic resonators with irregularly shaped boundaries 

. Even in the absence of pumping random media show 
interesting interference effects such as coherent backscat- 
tering or localization of light Q. The theoretical in- 
vestigation of the quantum optical properties of random 
media started with the demonstration of deviations 
from blackbody radiation in the case of one dimensional 
scattering. In a pioneering work |^ Beenakker gener- 
alized these results to multidimensional chaotic scatter- 
ing by establishing a general relationship between the 
photocount statistics of a linear random medium and its 
scattering matrix. Using that relationship, Beenakker 
and coworkers predicted an excess photon noise due to 
multiple scattering , computed the photon noise power 
spectrum and investigated the effect of photon lo- 
calization on the photocount statistics |ll|. 
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FIG. 1. Sketch of chaotic laser cavity with a partly trans- 
mitting mirror. The cavity is connected to a waveguide and 
a frequency selective photodetector. 



In contrast to the wealth of results for linear random 
amplifiers, little is known about the photon statistics of 
random media above the lasing threshold. This is an im- 
portant open problem since the experimental signatures 
of random lasing are still controversial Q| . In the present 
paper we address that problem for a chaotic laser res- 
onator in the regime of single-mode lasing. We show that 
the chaotic nature of the cavity modes gives rise to fluc- 
tuations of the photocount on top of the quantum opti- 
cal fluctuations known from laser theory. Chaos-induced 
fluctuations are found when a single-mode photodetec- 
tion is performed over an ensemble of modes. The en- 
semble may be obtained from a single resonator upon 
varying suitable parameters or from different resonators 
with small variations in shape. The factorial moments 
of the photocount display strong fluctuations from one 
mode to another. However, we show that the probability 
density for each factorial moment depends only on gen- 
eral symmetries of the ensemble and on four parameters 
describing the laser dynamics; these are the coefficients 
of linear gain, mean escape loss, absorption loss and sat- 
uration of the amplifying medium. 

A sketch of the system is shown in Fig. 1: A chaotic res- 
onator with an irregularly shaped boundary is homoge- 
neously filled with an active medium. A partly transmit- 
ting mirror reflects impinging radiation with mean prob- 
ability R back into the resonator. The fraction T — 1 — R 
of radiation is transmitted and injected into a waveguide 
supporting M transverse modes. A frequency selective 
photodetector counts the transmitted photons with effi- 
ciency 1. We consider the case of a quasi discrete cavity 
spectrum; then the mean cavity escape rate F = AITAuj 
is smaller than the spacing Alo of adjacent resonator 
modes. Accordingly, the mirror transmission coefficient 
T must be small, MT <C 1. The counting time t will be 
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assumed large enough that the radiation from individual 
cavity modes can be resolved, tAw 3> 1. 

As for the active medium, we allow for N active atoms. 
The /ith atom "sees" the field mode through a coupling 
constant which is proportional to both the atomic 
dipole matrix element and the value of the lasing mode at 
the location of the atom. We only consider the simplest 
situation in which the characteristic times for atomic 
pump and losses are short compared to the mean life 
time of a photon in the cavity. The effect of the atoms 
on the field mode can then be represented by three pa- 
rameters A, B, C that characterize the linear gain, the 
nonlinear saturation, and the total loss, respectively 
The first two of these depend on the coupling constants, 
A ^ 9% ^ gfi- Inasmuch as the resonator 

mode represents wave chaos, it varies quasi-randomly on 
the scale of the optical wavelength A Coupling 
constants associated with atoms separated by more 
than A therefore behave as independent random num- 
bers. The coefficients A and B, however, become sharp 
in the limit N ^ 1, due to the central limit theorem. The 
total loss rate C = F -|- k is the sum of the photon escape 
rate F and the absorption rate k. Here, the absorption 
rate may be considered fixed while the escape rate is a 
random quantity, as will be revealed presently. 

To calculate the photocount statistics we use the 
input-output theory of Gardiner and Collet |p^,p|. 
Waves entering and leaving the waveguide are described 
by Af -component annihilation and creation operators 
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a'"', a""*, a°"*^, which obey the commutation re- 



lations 

[ap{t),alit')] = 5p,Sit-t'), [ap{t),aqit')]=0. (1) 

Here, p,q = 1, . . . ,M label the transverse modes and 
a = a'" or a = a°"*. The boundary condition 



af{t)+a;^Ht)=j,b{t)., 



(2) 



connects incoming and outgoing radiation in each trans- 
verse mode with the annihilation operator b of the cavity 
mode. While the boundary condition (Q) results from as- 
suming a linear coupling between the waveguide and the 
cavity field, no restriction is imposed on the intracav- 
ity dynamics. Therefore, the input-output relation (|^) 
holds both below and above threshold. The outcoupling 
amplitudes "fp must be considered independent random 
quantities with Gaussian statistics, since they represent 
the local fluctuations of the resonator mode across the 
outcoupling mirror. In fact, their randomness constitutes 
the principle effect of the wave chaos within the resonator 
on the ensemble fluctuations of the laser output. 

In the linear regime below threshold, one can elimi- 
nate the cavity operators from Eq. and express the 
outgoing radiation in terms of the incoming radiation, 
the intracavity noise, and the S-matrix [[lUy]. One thus 



obtains the photocount fluctuations Q through the S- 
matrix statistics known for chaotic scattering. 

We proceed to the nonlinear regime near and above 
threshold. For the case of a vacuum input, the boundary 
condition allows to express all variances of the output 
field entirely in terms of the variances of the cavity field 
and the Gaussian statistics of the 7p. We characterize 
the photocount statistics through the factorial moments 
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m(m — 1) • • • (m — r + l)p(m) 



(3) 



of the probability p{m) that m photons are counted in a 
time interval t. That probability is given by 



p{m) = —I: W^e-^ :) 



•^0 p=l 



W 



(4) 



Here (• • •) denote the quantum steady-state average and 
the colons demand normal and a certain time ordering. 
Combining Eqs. with the input-output relation 

(B) for a vacuum input we get 



(5) 



where F = X^pii \lp\''^ the escape rate and the rth 
power of the integrated cavity field intensity 



dt'b\t')b{t') 



(6) 



Clearly, the first moment alias mean photocount, /ii = 
Tt{b%), is a purely static quantity as it is proportional 
to the stationary mean photon number inside the cavity. 

The factorial moments (||) specify the output statistics 
of a single cavity mode with the random escape rate F. 
The distribution of F over an ensemble of modes is the 
xi distribution 



P(F) = A^F''/2-iexp(-z/F/2F) 



(7) 



well-known from random-matrix theory ]15| . Here, v — 
(3M is an integer and A,y a normalization constant. The 
value of the parameter /3 depends on whether the sys- 
tem is time-reversal invariant {P — 1) or whether time- 
reversal invariance is broken (/3 = 2) The special 
case M — (3 — lis known as the Porter-Thomas distribu- 
tion. Together with F the factorial moments fj,r become 
random quantities. Their distribution is given by 



Vi^ir)^ / dTPiT)6{fir-r'-{:F:)). 



(8) 



Note that the right hand side involves a twofold average, 
the quantum optical average (represented by the brack- 
ets (•••)) and the ensemble average over the cavity modes 
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(represented by the integral with the probabiHty distri- 
bution P(r)). We emphasize that the quantum optical 
average (: :) depends on T through the total loss rate 
C. We now discuss the result (H) in various limiting cases. 

We first consider the case M ^ 1 of many transverse 
modes in the waveguide. The diameter of the waveguide 
is then much larger than the wavelength of the cavity 
mode. A simple saddle point argument shows that P(r) 
for large M approaches a Gaussian distribution with 
mean T ^ M and standard deviation AT = r/^/3M/2. 
The relative fluctuations are small, AT/T ~ 1/VM 
for M ^ oo, and the same is true for the fluctuations of 
all factorial moments. For cavities with large outcoupling 
mirrors we thus recover the sharp factorial moments one 
is used to from non-chaotic resonators. 

Second, we investigate the limit of vanishing photo- 
count, ^ 0- According to Eq. (||), this is the weak- 
coupling limit r ^ for which the photons in the cavity 
field can hardly escape into the waveguide. The laser dy- 
namics becomes independent of the outcoupling loss as 
the total cavity loss C is fully dominated by the absorp- 
tion loss. All cavity field moments become independent 
of r, and Eq. (||) reduces to fir ^ T^- Substitution into 
Eq. yields the power-law behavior 

Vifir) - fi'/'''^^^'^^-', (9) 

for ^ 0. A special case of this result is M = /3 = 1 
for which the distribution of the first and second moment 
diverge as and ^2 '^^^ fo'" Mi: — > 0, respectively. 

The third case is the short-time regime t ^ tc where tc 
is the correlation time of the intensity fluctuations of the 
cavity mode. During counting intervals that short the 
field intensity cannot vary appreciably. Therefore, the 
integrated intensity / becomes the product of t with the 
instantaneous photon number n = 6^6 in the cavity. The 
factorial moments of the output take the simple form 

fir^irtn-.n^:), (10) 

and one easily verifies that the normally ordered moment 
(: :) occuring here reduces to the rth factorial moment 
of the stationary photon number distribution of the lasing 
mode. The latter has the form [QJ|] 



Note that P„ depends on the three laser parameters A, 
B, and C; in particular, the nonlinearity B enters through 
the so-called saturation photon number ng = A/B. The 
symbol Af represents a normalization constant. 

The calculation of the mean photocount /ii reduces 
to a steady-state average for arbitrary counting time. 
We calculated fii numerically as a function of T using 
the distribution (|ll|). Substituting the result /ii(r) into 




0.2 0.4 0.6 0.8 1 



Ml 

fMnax 

FIG. 2. Distribution P(^i//imax) as a function of tlie di- 
mensionless mean photocount /ii//imax for /3Af — 1 and four 
sets of laser parameters. Rates are given in units of ^ = 1, 
the nonlineariti^is B = 0.005. (a) k = 0.7, T = 0.02 (solid 
line); k = 0.7, T = a 2 (dashed line), (b) k = 2.0, T = 0.5 
(solid line); k — 2.0, F — 4.0 (dashed line). The insets in (a) 
and (b) show as a function of T/A. The arrow indicates 
the threshold-value of F below which lasing takes place. 

Eq. (jsj) and carrying out the integration over F, we ob- 
tained Vipi). The results for time-reversal invariant cav- 
ities connected to a single-mode waveguide are plotted in 
Fig. 2, for f3M = 1 and four different sets of parameters. 
The two distributions in Fig. 2(a) correspond to .4 > C, 
i.e. to lasers above threshold in the ensemble average. 
By contrast, ^ < k for the distributions of Fig. 2(b); all 
lasers of those ensembles are below threshold irrespective 
of the escape rate F. We note that all distributions are 

— 1/2 

strongly non-Gaussian. They are all peaked as fj.-^ at 
small photocount, in accord with the general argument 
presented above for the asymptotics at ^ 0. Two fur- 
ther features spring to the eye and demand explanation: 
First, above threshold but not below we encounter an ad- 
ditional peak at maximum photocount; second, for cer- 
tain cases (dashed lines in Fig. 2) the distribution P(/ii) 
displays a shoulder for sub-maximal fii. 

The origin of these structures lies in the F-dependence 
of the mean photocount fii. That dependence is depicted 
in the insets and seen to be qualitatively different above 
and below threshold. While fj,i increases monotonically 
with F in the below-threshold case, it develops a max- 
imum at an intermediate value of F when the laser is 
above threshold. This behavior can be understood from 
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the simple analytic expressions 



C-A 



A-c 
c 



for C 
for C 



r + K<A 

r + K> A 



(12) 



that follow from the photon number distribution (|llh suf- 
ficiently far from threshold. According to Eq. (]l2[), /ii 
rises linearly with T out of the origin and approaches an 
asymptotic plateau for very large F, as visible in Fig. 2. 
The rise to the plateau is monotonic when A < k since 
the below-threshold case of ( 12 ) applies for all values of F; 
the maximum photocount /imax is then just the plateau 
value. The laser of Fig. 2(a) is above threshold for small 
F but below for large F. The maximum photocount //max 
then arises for an intermediate value F*. A simple argu- 
ment can be employed to determine the border between 
the below- and the above-threshold ensemble. The argu- 
ment follows from the observation that fii (F) approaches 
its plateau-value from above for the above-threshold en- 
semble while the plateau is approached from below when 
the ensemble is below threshold. From the lower case of 
Eq. (p^, which becomes exact for F ^ oo, this yields 
the threshold condition A = k. To estimate the value of 
F* we may employ the above-threshold case of Eq. ( p^ ) 
and find F* = V An — k. 

Based on this understanding of /ii(F) one can appre- 
ciate the above-threshold peak of 'P(/ii) at /ii — /imax- 



Mmax + 5M"(r*)(F 



F*)^ 



the 



Substituting /ii(F) 
photocount distribution (|^) we find that V{^i) has a 
square-root singularity |/ti — /fmaxl^^^^ which is precisely 
the peak depicted in Fig. 2(a). Clearly, no such peak can 
arise in the below-threshold case of Fig. 2(b) as the pho- 
tocount increases monotonically with F. Inspecting Eqs. 
(0); ® the limit F ^ oo we rather find that P(/ii) 
vanishes when /ti approaches /imax- 

The shoulders at sub-maximal photocount are caused 
by amplified spontaneous emission below the laser thresh- 
old. Formally, the shoulders arise from the asymptotic 
plateau of /ti for large F. The definition (^) immedi- 
ately implies enhanced probability for photocounts /ii 
in the vicinity of that plateau. Note that the shoulder 
is invisible for the distribution shown as a solid line in 
Fig. 2(a), since for this distribution large values of F are 
strongly suppressed by the small mean value F. Further 
note that the shoulders for both curves in Fig. 2(b) lie 
closer to /imax than for the dashed curve in Fig. 2(a) 
since /t — > /tmax coincides with F — > oo in the regime 
below threshold. 

In contrast to the mean photocount which can be ex- 
pressed in terms of the stationary distribution of the 
laser, all with r > 2 involve dynamical information 
through correlation functions with r — 1 time arguments. 
We defer the discussion of the higher factorial moments 
and their distribution to a separate publication p^ . 

We finally compare our results with related fluctuation 
phenomena in other areas of physics and discuss possi- 



ble experimental tests. In nuclear physics the Porther- 
Thomas distribution describes level-width fluctuations 
in neutron scattering ||l^. The amplitude fluctuations 
of Coulomb blockade oscillations in semiconductor quan- 
tum dots are also of the Porther-Thomas type |l^. In 
both cases, as well as for the chaotic lasers studied in this 
paper, the fluctuations result from the chaotic nature of 
wave functions. However, in chaotic lasers new interest- 
ing features arise due to the interplay of wave chaos with 
the nonlinear dynamics of the laser. As a consequence of 
this interplay, the distribution of the mean photocount 
can strongly deviate from the Porter-Thomas distribu- 
tion. To test the predicted mode-to-mode fluctuations 
experimentally, one must study the photocount statistics 
for an ensemble of chaotic modes. It seems feasible to 
generate such ensembles in tunable lasers e.g. by shape 
variations, or by the injection and displacement of artifl- 
cial scatterers in the case of microwave cavities. 
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